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ABSTRACT 



Multiwavelength observations of blazars such as Mrk 421 and Mrk 501 show 



o 
o 
o 

fS| ■ that they exhibit strong short time variabihties in flare-hke phenomena. Based on 

• the homogeneous synchrotron self-Compton (SSC) model and assuming that time 

. variability of the emission is initiated by changes in the injection of nonthermal 



electrons, we perform detailed temporal and spectral studies of a purely cooling plasma 
system, using parameters appropriate to blazars. One important parameter is the total 
. ^ injected energy £ and we show how the synchrotron and Compton components respond 

' as £ varies. When the synchrotron and SSC components have comparable peak fluxes. 



we find that the SSC process contributes strongly to the electron cooling and the whole 



fSj ■ system is nonlinear, thus simultaneously solving electron and photon kinetic equations 

. is necessary. In the limit of the injection-dominated situation when the cooling 

. timescale is long, we find a unique set of model parameters that are fully constrained by 

I observable quantities. In the limit of cooling-dominated situation, TeV emissions arise 

Qh| mostly from a cooled electron distribution and Compton scattering process is always 

O ' in the Klein-Nishina regime, which makes the TeV spectrum having a large curvature. 

' Furthermore, even in a single injection event, the multiwavelength light-curves do 

■ not necessarily track each other because the electrons that are responsible for those 

emissions might have quite different lifetimes. We discuss in detail how one could infer 
^ , important physical parameters using the observed spectra. In particular, we could infer 

the size of the emission region by looking for exponential decay in the light curves. We 
could also test the basic assumption of SSC by measuring the difference in the rate of 
peak energy changes of synchrotron and SSC peaks. We also show that the trajectory 
in the photon-index-flux plane evolves clockwise or counter-clockwise depending on 
the value of £ and observed energy bands. 
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1. INTRODUCTION 

Blazars are a class of flat radio spectrum, core-dominated active galactic nuclei (AGNs). The 
overall radiation spectra of blazars show two broad peaks in the vFi, space; one is between IR and 
X-rays, and the other in the 7-ray regime (e.g., von Montigny et al. 1995|) . Flares also have been 



observed at X- and gamma-ray bands by multiwavelength observations of Mrk 421 (e.g., Macomb 
et al. 1995; Macomb et al. 1996 for erratum; Buckley et al. 1996) and Mrk 501 (Catanese et al. 
1997; Plan et al. 1998). The tremendous luminosity and fast time variabilities from blazars have 
led to the usual arguments that relativistic motion is occurring in the emitting plasma. Moreover, 
the favored scenario to explain these sources is that we are viewing nearly along the axis of a 
relativistically outflowing plasma jet that has been ejected from an accreting super massive black 
hole (e.g., [Blandford fc Rees 19781). 



Although the origin of these multiwavelength spectra is still under debate, several models 
on the radiative processes have been put forth, in particular, models of Compton scattering of 
synchrotron photons or external photons have been developed in recent years (e.g.. Bloom & 
Marscher 1996; Inoue & Takahara 1996; Ghisellini & Madau 1996; Dermer, Sturner, & Schlickeiser 
1997; Mastichiadis & Kirk 1997; Sikora et al. 1997; Bottcher, Mause, & Schlickeiser 1997; 
Georganopoulos & Marscher 1998; Ghisellini et al. 1998). Most of these calculations are either 
semi-analytic, or for steady state situations, or not including the Compton scattering process 
self-consistently. The main purpose of this paper is to improve upon this situation. The physics 
of how energy is dissipated into relativistic particles is, unfortunately, not well understood (see. 



however, Romanova &: Lovelace 1997) and will not be treated fully in this paper. 



Among various blazar models, synchrotron self-Compton (SSC) models have received a fair 
amount of attention, by virtue of its simplicity and its possible predictive power. In these models, 
it is proposed that the nonthermal synchrotron emission forms the radio-through-X-ray continuum, 
and that the Compton scattering of these (soft) synchrotron photons by the same nonthermal 
electrons produces the gamma rays (~ GeV - TeVs). In this paper, we focus on the so-called 
homogeneous SSC model where a spherical blob of uniform relativistic plasma is postulated. 
Even with such a greatly simplified picture, a number of parameters have to be invoked, whose 
interplay gives rise to a rich dynamic behavior of the observed radiation. Of particular interest is 
the correlated variabilities in X-ray and 7-ray fluxes, since they represent the tail of nonthermal 
electrons which have the shortest cooling timescale. Although the generic multiwavelength spectra 
from radio to TeV can be fitted by a steady state model with fixed parameters (e.g, Kataoka et 
al. 1999), time-dependent calculations almost always offer stronger constraints. Furthermore, 
when the self-Compton component contains a comparable or even larger fraction of the radiative 
energy than the synchrotron component, the whole problem becomes inherently nonlinear and 
both components need to be calculated simultaneously and self-consistently. This naturally leads 
to the need of solving coupled, time-dependent, nonlinear particle and photon kinetic-equations. 
Moreover, by examining the energy-dependence of flare data at gamma-ray energies, one could 
potentially discriminate between SSC and external Compton-scattering origins of the seed photons 
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(Dermer 1998). 

The simplest model for time variability of blazars (Mastichiadis & Kirk 1997; hereafter 
MK97) assumes that electrons obeying a power-law distribution are injected uniformly throughout 
a relativistically moving blob over an extended period of time, and that electrons cool by 
both synchrotron radiation and Compton scattering. The blob is assumed not to accelerate or 
decelerate, and the energy loss by Compton scattering of photons impinging from outside the blob 
is assumed to be small in comparison with the synchrotron self-Compton loss. MK97 reproduced 
the qualitative behavior of the energy-dependent lags and the hysteresis diagrams (Takahashi et 
al. 1996). Much of the work presented here follows closely to the previous study by MK97, but 
we are using a completely different kinetic code which will be discussed in later sections. Kirk, 
Rieger, &: Mastichiadis (1998) further modeled the evolution of synchrotron emission, calculating 
acceleration and cooling regions separately, though Compton scattering was not included. 

In this paper we present a detailed study of the time-evolution of an electron-photon plasma 
(the positive particles could be either protons or positrons) by solving the kinetic equations 
numerically. We briefly describe our model in §^ and show numerical results in §^. Summary is 
given in §Q. 



2. MODEL 

We assume that observed photons are emitted from a blob moving relativistically towards us 
with a Doppler factor V = [r(l — /3r cos 0)]~^, where F is the Lorentz factor of the blob, /5rc is the 
velocity with c being the light speed, and 9 is the angle between the direction of blob motion and 
the line of sight. The blob is a spherical and uniform cloud with radius R. Relativistic electrons 
are injected into the blob and produce high energy emission. Electrons and photons are uniformly 
distributed throughout the blob. The most important physical processes for electrons include 
synchrotron radiation and Compton scattering. The spectra of electrons and photons in the blob 



are calculated by solving the kinetic equations described below (see also Coppi fc Blandford 1990 ). 



The kinetic equation describing the time-evolution of electron distribution n(7) is given by 



dn _ d 

dt d"f 
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where n is the electron number density per 7, 7 is the electron Lorentz factor, and tg.esc is the 
time for electrons to escape from the blob. The term (d7/dt)ioss represents various electron energy 
loss processes, such as synchrotron and Compton scattering; their corresponding energy diffusion 
is given as Dg. We also include the synchrotron "boiler" effect (Ghisellini et al. 1988) and 
other processes such as Coulomb collisions, though they are not important in the present study. 
Pair-production and annihilation are not treated in the present code, though they tend to be less 
important too. In our code, the particle equation is actually discretized in the momentum space 
so that thermal particles and their processes can be handled accurately. This is less important for 
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AGN jet parameters but will be very useful for modeling the emissions from stellar-mass black hole 
systems. Note that equation (|^) assumes a continuous electron energy loss (i.e., Fokker-Planck 
type). This assumption does not account for the situation when there is a significant energy loss in 
a single Compton scattering which is important in the Klein-Nishina regime. This discrete nature 
of the Compton energy loss is, however, included in the photon kinetic equation, i.e., equations (|2|) 
and (^). We have checked the accuracy of the continuous approximation of the Compton energy 
loss in equation (1) in view of energy conservation. Our numerical tests show that the total energy 
is conserved with the accuracy better than 5 per cent after 10 R/c, when the SSC component is 
dominant and the scatterings occur frequently in the Klein-Nishina regime. 

The relevant kinetic equation for the time-evolution of photons is given by 

— ^- — = nc(e) + nem(e) - n^hs{e) - -f > (2) 

C'' '•ph,csc 

where nph(e) is the photon number spectrum per unit volume per unit energy e. Compton 
scattering is calculated as 

nc(e) = -nph(e) J d-f n{-f) Rc{e,^) + J j de' P{c,e' Rc{e\-i)n^^{e')n{^) . (3) 

First term of the right hand side of equation (^) denotes the rate that photons with energy e 
are scattered by electrons with number spectrum 71(7) per unit volume per unit 7; Rc is the 
angle-averaged scattering rate. Second term of equation denotes the spectrum of scattered 
photons: P(e;e',7) is the probability that a photon with energy e' is scattered off by an electron 
with 7 to have energy e. The probability P is normalized such that / P{e\e' de = 1. The 
details of Rc and P are given in Jones (1968)| and the appendix A of Coppi &: Blandford (1990). 



We use the exact Klein-Nishina cross section in the calculations of Compton scattering. Photon 
production and self-absorption by synchrotron radiation are included in nem(e) and ?iabs(e)) 
respectively. The synchrotron emissivity and the absorption coefficient are calculated using the 
approximations given in Robinson and Melrose (1984) for transrelativistic electrons and Crusius 
and Schlickeiser (1986) for relativistic electrons. External photon sources such as disk photons 
can be included, though they are not considered here. The rate of photon escape is estimated 
as ?T-ph(e)/iph,esc- Since we are in the optically thin limit, we set tph,esc = R/c, which is a good 
approximation. The photon spectra from solving equation (^) has been extensively compared 
with those from Monte-Carlo simulations and we have found very good agreement between them 
(Kusunose, Coppi, & Li 1999). 

The comoving quantities are transformed back into the observer's frame using eobs = e ^/ (l+'Z) 
and dtobs = dt{\ + z)/T>, where z is the redshift of the source. 

We assume that electrons are injected obeying a power law in energy: 

Q(7,t)=Qo(t)7-^ 

with 7min < 7 < 7max- The total energetics of the electrons can be represented by a compactness 
parameter which is proportional to L/R, where L is the source luminosity. We do not consider 
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specific acceleration mechanisms in this paper. Thus particles are just being injected into the 
blob over a finite time. We emphasize that in order to be consistent with the basic assumption 
of spatial homogeneity, we require the injection time tinj to be longer than tdyn = R/c. In effect, 
we can not probe variabilities shorter than tdyn in the comoving frame. Other physical effects 
such as adiabatic loss via expansion might play an important role but is not considered here. It 
is unfortunate that we need such a large number of parameters to proceed with the calculations, 
and this is the main reason we opt not to add further complications such as acceleration. 

3. RESULTS 

The dynamic behavior of the emission spectra is controlled by several timescales, namely, 
the cooling time tcoob the dynamic time tdym the injection duration finj, and the escape time 
icsc- The causality argument requires that t^yn < tinj,icsc) whereas tcool can be smaller than tdyn- 
For SSC models, both synchrotron and Compton processes contribute to the electron cooling, so 

l/*cool = 1/isyn + l/^ssc, where tsyn = 7/l7syn| and issc = T/ITsscI- 

In the following analyses, we will divide our results into two major parts, based on whether 
tcooi is longer or shorter than t\^y In the limit of tcooi ^ *inj > tdyn-, where icooi is evaluated 
using the highest particle energy, the injected particle distribution does not change appreciably 
during the injection process. We call this the injection-dominated limit. On the other hand, if 
t'm] > *dyn > ^cooh then particlcs will be sufficiently cooled while the injection still occurs, and 
emissions are from a cooled particle distribution rather than from the injected one. We call this 
the cooling-dominated limit. Consequently, we expect rapid variations in both fluxes and spectra 
in the latter case and relatively slow spectral variations in the former case. 

The primary purpose of this paper is to understand the dynamics of SSC model. We thus 
have chosen a broad range of parameters rather than try to fit any specific source spectrum, but 
we certainly use parameters thought to be applicable to those sources (Mrk 421 in particular) to 
guide our calculations. 

3.1. Long Cooling Time Limit 

We show that, in this limit, a set of SSC model parameters can be uniquely determined from 
the observable quantities. We use observations of Mrk 421 as an example and further discuss their 
implications for multiwavelength observations. 
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3.1.1. A Unique Solution 

Let £ represent the energy (in ergs) injected in nonthermal electrons, which we assume can be 
described as A^e(7) = 7"^ and < 7 < 7rnax. Then iVo = (2 - p)£ /[meC^{-('^^l - 7^7n )] (for 
p ^ 2). In the limit of long cooling time, we can use the injected electron distribution to calculate 
the synchrotron and SSC fluxes. Thus, for the peak energies of synchrotron and SSC, we have 

2?7Lx^ « X = , (4) 

^?7max ~ l^ssc/1.236 X 10^0 = i^^sc , (5) 

where Vsyn and u^^c are the synchrotron and 7-ray peaks in Hz, respectively, and B is magnetic 
fields in Gauss. The numerical normalization factors are easily obtainable using the standard 
expressions for the peak synchrotron energy and inverse Compton peak energy in the KN limit. 
For Mrk 421, we have v^yn ~ 1-72 x 10^^ (~ 2 keV) and i>ssc ~ lO'^ (~ 5 TeV). From equations (|) 
and (III), we get 

B = VD,^jDl^ . (6) 

The relative ratio r] of SSC to synchrotron fluxes for Mrk 421 is close to 1. This ratio can be 
approximately represented by the ratio of the comoving-frame synchrotron photon energy-density, 
J/syn = Lsy„/(47ri22cp4)^ to the magnetic field energy-density, Ub = B'^/{8tt). But the exact value 
of 1] could be quite different from this estimate owing to several factors: the end point effect at the 
peak of both synchrotron and SSC fluxes, |^ the KN effect of Compton scattering for producing 
TeV emission, and the fact that the electron distribution will be slightly cooled even though the 
cooling timescale is long. It is very difficult to get an exact analytic value to account for all these 
effects, so we introduce a correction factor fc in calculating 77. Thus, we have 

V = fc^ or Lsyn « VUttR^c Ub {^/fc) • (7) 
Ub 

Furthermore, since all blazar sources are observed to be highly variable, an additional constraint 
has usually been proposed that 

R^VcU^, , (8) 

where tvar is defined as a variability timescale in the observer's frame. 

Combining all the equations given above, we have 4 equations for 4 independent variables 
(i.e., P, i?,7max7 and R). The supplemental information include z^syn, t'ssc, -^syn, and tvar, with a 



^For p < 3, the peak of i^L^ is at Usyn oc 7^3x^1 but its flux is smaller than that obtained using the 5— function 
approximation to the scattering cross section. To be consistent with our numerical code results which have used 
the exact synchrotron emissivity formulations, we have introduced a reduction factor of /syn(< 1) in our simplified 
analytic estimates for the peak synchrotron fiux (see also equation (jl^). Note that this reduction only applies to the 
end points of emissivity. For energies much smaller than Vsyn, the 5— function approximation is quite accurate. 



-7- 



somewhat variable factor fc- Thus we can uniquely determine a solution set for all the parameters. 
The Doppler factor from this solution (denoted by a subscript 's') can be expressed as 



The peak luminosity (Lsyn) of the synchrotron component in the observer's frame can be 
estimated using vLy at z/gyn (again taking into account the end-point effects), 

^syn ~ /syn^^^^-^0-TC?7B— ^ 7max , (10) 

where o"t is the Thomson cross section, /syn(< 1) represents the reduction of synchrotron flux at 
the end point, is due to the Doppler boosting. Observationally, Lgyn at ~ 2 keV is > 6 x 10^^ 
ergs s~^ with a luminosity distance of 3.8 x 10^^ cm {z = 0.0308) for a = 1/2 and A = 
cosmology. Here, Hq = 75 km s~^ Mpc~^ is assumed. Also, fitting of Mrk 421's synchrotron 
spectrum suggests p ~ 1.65 (MK97). From this, we can determine the rest of the parameters for 
this unique solution (again denoted by a subscript 's') 

Bs = P.Z>syn/^L ~ 0.038 G , (11) 

7max,« = V^sc/Vs « 4.5 X 10^ , (12) 

= Ctvar^s ~ 6.5 X 10^^ cm , (13) 

1.55 X 10^ -Z^syn . 1 ,„46 Iiax 

^« = (2-»)f 7 ^^4.1x10 6 ergs , (14) 

\^ l-'JJsyn /max, s-'-'s -^s 

where we have used p = 1.65 and /gyn = 0.4. Q More importantly, we can check our original 
assumption that cooling time is long compared to tdyn = R/c. This is obviously satisfied since 
isyn ~ (6vrmec/(TT)/(7max-B^) « 1.2 X 10^ sec, which is much longer than tdyn ~ 2.2 x 10^ sec. 

Using the above derived parameters, we solve equations (|^) and (|2|) simultaneously and follow 
the evolution until 20tdyn- A total energy of £ is injected in nonthermal electrons over a comoving 
timescale of tjnj = 2tdyn- In these calculations electrons are not allowed to escape. Figure |l] shows 
the time evolution of electron and photon distributions as the system evolves. Note that the time 
for synchrotron and SSC components to reach their peak fluxes is different, and that it happens 
after the electron injection has stopped. To qualitatively understand this, we can write the photon 
kinetic equation symbolically as 

— = Production(e) — Escape(e) . (15) 



^The fact that fc and fsyn are both chosen as 0.4 is a coincidence. The evaluation of fc involves end-point effects 
from both synchrotron and Compton scattering, whereas /syn is only concerned with synchrotron. The value 0.4 is 
obtained by comparing the analytic estimates with the exact numerical calculations. 
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Thus, the photon flux at energy e will increase if the production rate is larger than the escape and 
decrease if escape is quicker. This determines when the peak of photon flux at certain energy is 
reached. Using photon flux at keV as an example, the production of these photons still continues 
even when the electron injection stops because the cooling timescale is much longer than R/c. 
Eventually as electrons cool, they can no longer produce keV synchrotron emissions, the flux at 
keV starts to decline. 

Figure |2| shows the light curves at different energy bands expected from this injection event. 
Since the cooling time is rather long, only the energy bands corresponding to the tail of electron 
distributions show short time variabilities due to that electron injection is turned on and off; 
whereas other energy bands show a long plateau, representing a balance between the photon 
production and escape. A clear prediction from this is that there should be very little spectral 
evolution except at the peaks of synchrotron and SSC. All these are commensurate with the 
dynamics of electron cooling. 

In order to further differentiate the role of synchrotron versus Compton cooling on electrons, 
we plot the ratio of |7ssc/7syn| as a function of electron energies at different times in Figure ^. 
It is clear that SSC cooling becomes more important than synchrotron cooling when the photon 
energy density builds up within the system as time proceeds. After reaching the peak, the SSC 
cooling starts to decrease as the photon energy density decreases. The dependence of this ratio on 
the electron energy is partly due to the KN effect. This figure clearly indicates that one can not 
ignore SSC cooling in estimating certain parameters and that the evolution is very nonlinear, thus 
a self-consistent numerical calculation is required for fitting the data more accurately. 

Note that the spectra given in Figure || is not intended to be an accurate fit to the observed 
spectra from Mrk 421. In fact, the TeV spectrum of Mrk 421 is known to be roughly a power law 
( Krennrich et al. 1999| ), but the generic spectrum in TeV obtained here has a clear curvature, due 
to the fact that the Compton scattering is in the KN regime. Nevertheless, this exercise allows us 
to establish a parameter space where a reasonable fit to the actual spectrum might be obtainable, 
and it has the nice feature that the electron energy distribution retains the injected form without 
much softening, which greatly simplifies the analysis. 



3.1.2. Parameter Variations on the Unique Solution 

In this subsection, we explore how sensitive the above results are to parameter variations. 
The parameter we want to emphasize is the ratio r/ of the SSC component to the synchrotron 
component. From equations (^) and (P), and holding other parameters unchanged, we can see that 

r? oc V-^ oc , (16) 

which implies that a small change in magnetic field and/or Doppler factor could result in a 
large variation in r/. Figure ^ shows this effect when B is being varied. The general trend from 
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equation (16) is indeed confirmed, i.e., smaller/larger B (versus Bg) gives much larger /smaller rj. 
For B < Bg, the variation amplitude in rj does not exactly follow equation (p!6|). We attribute 
this discrepancy mostly to variation in because both 7max and synchrotron photon energy 
in comoving frame are varying for different B. Additionally, when r/ > 1 (SSC cooling is more 
important than synchrotron cooling), electron distribution is subject to stronger cooling. This is 
evident in its rapid spectral variation in the middle panel of Figure 0. For larger B which results 



in ?7 < 1, equation (16) is mostly confirmed. Note that the Doppler factor P, using equation 



is now getting uncomfortably large. 

To conclude, in the limit that the cooling time of the highest energy electrons is longer than 
the dynamic timescale on which injection occurs, we can most likely find a unique set of parameters 
that roughly satisfy the observational constraints. The most important prediction is that even 
though the fluxes of the synchrotron and SSC peaks can vary by a large factor (> 10) in a short 
time whose time histories look like "pulses", the duration of emission at other wave bands (such 
as GeV, MeV, eV) can be considerably longer by at least a factor of 4, which is commensurate 
with long electron cooling time at those energies. 

This parameter space has some interesting implications for interpreting multiwavelength data 
from blazar monitoring campaigns. Using keV and eV bands as an example (similar arguments 
can be made for TeV and GeV bands too), since the lifetimes for keV-producing and eV-producing 
electrons are different, there is really no reason to expect their light curves to track each other 
and to have the same rise and fall patterns or timescales. Furthermore, when there are multiple 
injections occurring over a timescale shorter than the lifetime of keV-producing electrons but 
longer than the lifetime of eV-producing electrons, fluxes at keV could vary rapidly to reflect the 
multiple injections. Fluxes at eV, however, might only show a continuous increase with no obvious 
decline because of the accumulation of eV-producing electrons from multiple injections (see also 



3.2. Short Cooling Time Limit 

In this subsection we explore the limit where tcooi(7max) < ^dyn, i-e., electrons are efficiently 
being cooled while they are injected into the system. All the runs in this section have size 

= 1.5 X 10^^ cm which gives a comoving dynamic timescale R/c = 5 x 10^ sec, particle's 
7min = 10, 7max = 10^, and index p = 2. The Doppler factor is chosen as 10 and magnetic field 

= 0.1 G, though most of the conclusions depend weakly on D and B in this section. The 
particularly attractive feature of this limit is that it is possible to achieve short time variabilities 
for many wave bands, contrary to the case shown in Figure ^ from the previous section. 

One serious problem in comparing the theoretical results with the actual observations is that 
most observations need to accumulate over certain time interval (to collect enough photons) and 
different integration times are needed for different energy bands. So, in order to make a direct 
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comparison with the observations, properly averaged fluxes are needed, rather than the prompt 
flux we have presented above. This inevitably introduces many more additional parameters in 
determining how photons at different energy bands are sampled. To avoid these complications, 
as before, we will continue to present the prompt flux results, leaving the problem of integrated 
fluxes to future studies on detailed spectral fitting of particular sources. 

To quantify the relative importance of SSC versus synchrotron, the previous expression 
for rj has to be modified. To recapitulate, the comoving synchrotron photon energy-density is 

TT ~ 



[meC^ / (Att R'^ c)] Z^™'"' A''e(7)|7| dj, where the electron energy loss-rate through synchrotron 



radiation is given by 7 = —[Aa^ / {3meC)]UB^ ■ Thus, we get 



V = fc 
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Ub 
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(17) 



where the last expression applies to the case p = 2. On the other hand, if idyn > icooh a cooled 
electron distribution has to be used when calculating the synchrotron photon energy-density. A 
very rough estimate of the electron break energy 7br, beyond which cooling dominates can be 
given as 

(\7rm 1 

max (tinj , tdyn ) • 



CTT 7br-B^ 

The cooled electron distribution has the origina,! ind.Gx — j> between ^min 

and 7br, and roughly 

—•p — 1 between 7br and 7max- Since the number density of electron at high energy end is very 
small for p > 1, the previous expression for iVo (§ P^ ) still apphes. Then we can derive another 
expression for r] as 



?7c = fc 



37ri?2 



■P 



3-p 
7br 

/max 



+ 1 



fc 



37ri?2 



7br + ln(7max/7br) 



ln(7max/7mm) 



(18) 



where again, the last expression applies to the case p = 2. Note that this ratio rjc depends on 



B through 7br. Using the same parameters given previously, we find that equations (17) and 
(|T^ give r]c ~ (7br/7max)^- In other words, in order to reach the same relative ratio between 
synchrotron and SSC, more energy is needed (by a factor of 7max/7br for p = 2) if electrons are 
cooled efficiently during injection. 



3.2.1. Dynamics of a Single Injection 

In this subsection we concentrate on the dynamics of a single injection event lasting tjnj and 
its corresponding evolution of electron and photon distributions. The idea is to mimic individual 
flaring events in blazars and gain some basic knowledge of how synchrotron and SSC components 
are dynamically linked. To simplify the calculations and analysis, we choose 6 total injection 
energies with a factor of 10 increase from 10"^^ ergs to 10^^ ergs. These energies in nonthermal 
electrons are injected over a comoving timescale of tinj = 2tdyn- We solve equations (|l|) and (§) 
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simultaneously and follow the evolution until lOtdyn- Electrons are not allowed to escape. With 
these parameters, tcooi(7max) is much smaller than tdynj so electrons are appreciably cooled during 
injection. 

Figure |5| shows the i^Fi, spectra of these 6 injections taken at t = tjnj and with fluxes divided 
by S/W^"^. It is evident that synchrotron is the dominant cooling process for £ < 10^"^ cases, 
whereas SSC becomes the dominant cooling process as the photon energy density builds up in the 
10"^^ and 10"^^ cases. In fact, electron cooling is so significant in the 10 case that the maximum 
synchrotron flux is reduced by almost a factor of 10 compared to other lower injection energy 
cases, and its SSC peak energy is also much softer than others. 

As given in equation ([T^), the ratio of SSC to synchrotron is roughly proportional to 6. Thus 
as £ increases, so is this ratio. This is shown as the (almost) linear increase in the SSC peaks of 
Figure |5|. So we conclude that so long as the peak in SSC is less than synchrotron, the magnitude 
of increase in SSC peak will be the square of the magnitude of increase in synchrotron peak. When 
the SSC component becomes comparable to synchrotron component, the system becomes highly 
nonlinear, the estimate of 7br based on pure synchrotron cooling is no longer valid, even though 
equation (|l8|) probably still applies as long as a cooled electron distribution is used. 

A further point regarding the relative ratio of SSC versus synchrotron components is that the 
initial ~ 10 GeV - TeV production via SSC is in the KN regime, which reduces the SSC flux. This 
implies that an even larger £ is needed than those given in equation (^). Observations of Mrk 
421 and 501 seem to indicate roughly the same heights of synchrotron and SSC peaks. Thus, in 
modeling the time-dependent (or even steady state) emissions from these objects, a full KN cross 
section has to be used, as was done here in our code. 



3.2.2. Dynamics and Light Curves 

We now study in detail the full time evolution of three injection cases: £ = 10'^^, 5 x 10^'^, and 
-|^q49 gj.gg_ Xhey are shown in Figures ^ 0, and |8|, where the time-evolution of electron distributions 
and photon spectra are presented. Figures |9|, |l^, and 11 show the corresponding light curves 
of different photon energies for the above three injection cases. To qualitatively understand the 
light curves, we refer to equation (^) again. The peak of the light curves is reached when the 
production and escape are balanced, which depends on whether particle distribution has softened 
enough. Furthermore, once the production at certain photon energy has stopped, the pure escape 
process will produce an exponential decay. Since the photon escape timescale is tdyn = R/c, one 
could get an estimate of the size by fitting the decline portion of the light curves. This can be 



done using Figures y, |10|, and [11| where fluxes have been plotted in logarithm. The straight lines 
give a clear representation of photon escape, especially when it shows up in several energy bands, 
thus this might be a useful method in analyzing the real blazar data. 
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3.2.3. Spectral Variations 

The flux changes depicted in the above subsection are accompanied by large spectral variations 
too as shown in Figures |^ ^, and ^. These curves contain a wealth of information, some of which 
are rather parameter dependent. Nevertheless, we can draw some general conclusions: 

(1) Since the cooling time at 7max is the shortest timescale in our system, the electron 
distribution at high energy end is always substantially softened. Furthermore, the production of 
photons > tens of GeV is always in the KN regime using SSC model. Both effects make the TeV 
spectrum very soft, with a large curvature. [Additional effects such as intrinsic absorption at the 
source or intergalactic absorption by infrared background can cause further curvature (e.g., Coppi 
&L Aharonian 1999).] This curvature is not consistent with the TeV spectra we have seen from 
Mrk 421 ( [Krennrich et al. 199S| ), but it is perhaps consistent with the observations of Mrk 501 
where a curvature in the Compton component is clearly seen (Djannati-Atai et al. 1999). 



(2) The "hysteresis" in the relation of photon energy flux and spectral index was flrst pointed 
out by Takahashi et al. (1996) at the keV band. In Figures 13, and 14, we show the evolution 
of photon index as a function of energy flux in the observer's frame. Clockwise rotation is always 
seen at 1 GeV, regardless of whether synchrotron or SSC cooling dominates. Clockwise rotation 
is found at 2 keV for the case where synchrotron losses dominate the electron cooling {£ = 10^^ 
ergs). We flnd that this clockwise rotation at 2 keV is true for all cases with injection energy less 
than 10^® ergs (with the above injection form). This is mostly related to the fact that if we can 
associate 2 keV with the synchrotron peak, the spectrum always softens when its flux is decreasing 
because the electrons that can produce 2 keV photons have diminished. When the injection energy 
is large, the SSC loss is dominant {£ = 10^^ ergs), the hysteresis diagrams rotate in the opposite 
sense at 2 keV. The hardening at later time is actually due to the fact that 2 keV flux is now 
from the flrst generation of SSC, not synchrotron anymore. Different behaviors of the hysteresis, 
however, are found at 100 keV: counter-clockwise for £ = 10'^^ ergs and clockwise for £ = 10^^ 
ergs. Given these large variations and their sensitivity of various parameters, it is difflcult to use 
these "hysteresis" diagrams to draw flrm conclusions. 

(3) If observations show a synchrotron peak in 1 - 10 keV band, one should be very careful 
with fltting the spectrum in the 100 keV energy band (such as OSSE and BeppoSAX), since it 
is right in the region where synchrotron and SSC meet. There is a very large and rapid spectral 
evolution during the flare (e.g.. Figure p^). 

(4) In the rising part of the vFy spectra, such as 1 - 10 eV and MeV - GeV, the spectral 
index variation is much slower than the keV and TeV energy bands even though their fluxes vary 
by a large factor (e.g., Figures |l^, ^ and p^). 

(5) To further quantify the spectral evolution, we plot the peak energies of synchrotron and 
SSC components in vFy as a function of time in Figure ^ for the case with £" = 5 x 10^^ ergs. An 
important part is the early softening stage, where the synchrotron peak energy decreases as 7^ 
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whereas the SSC peak energy goes down first as 7 because the scattermg is in the KN regime. 
This effect might be observable with the current keV and TeV observational capabilities. Such a 
"correlated" evolution in peak energies might provide a definitive test of SSC. 

(6) As shown in Figure ^, when the SSC cooling becomes comparable to or dominant over the 
synchrotron cooling, the synchrotron peak becomes broader than those dominated by synchrotron 
cooling only. 

Some of the above conclusions might be testable using the current data collected on blazars, 
and some might require much higher quality data. 



3.2.4- Time- Dependent Injection 

In this subsection we show how different injection profiles change the light curves in a 
single injection event. Different from the previous subsection where a constant injection is 
used (box-shaped injection), we calculate another case with a linearly increasing injection rate 
(triangle-shaped injection), i.e., Q{'y) = Qo7^^V^dyn for < t < 2t(jyn, and (^(7) = for 
t > 2idyn- The time evolution is followed until t = lOtdyn and electrons are allowed to escape with 
te,esc = Stdyn- We use the same number of particles and same amount of total injected energy in 
both box- and triangle-shaped injections; the injected energy is 5 x 10'^'' ergs per 2tdyn in the blob 



frame. In Figure 16, we compare their light curves at 1 keV. The light curve for the box-shaped 
injection is asymmetric, because more electrons are injected at an early stage than in the triangle 
injection. On the other hand, the light curve from triangle injection shape is almost symmetric. 
In both cases, light curves decay exponentially after the end of the injection. 



3.2.5. Multiple Injections 

We now move to study other injection profiles, which are done by artificially turning the 
electron source term on and off. This is admittedly quite artificial. The purpose is to understand 
whether there are any generic features associated with these multiple injections, which might aid 
us on modeling the multiple flares often observed from blazars. In all the following runs, we allow 
electrons to escape with te,esc = 5tdyn- Other parameters are the same as the single injection case. 

Successive flares can be produced by repeated injections of nonthermal electrons. As an 
example of this picture, we present light curves for two repeated injection cases of nonthermal 



electrons. The top panels of Figures 17 and |18| show the injection profiles, which consist of two 
triangle injections separated by a long (Stdyn) and a short (2tdyn) intervals (in the comoving frame), 
respectively. In both cases 5 x 10"^^ ergs are injected in each "flare" with a triangle-shaped time 



profile. The lower panels of Figures 17 and 18 show the expected light curves calculated by our 



code. The shape of the light curves from each injection is very similar to that in Figure ^ where 
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a single injection is involved (i.e., quasi-sjnumetric). The peak fluxes of light curves in multiple 

injections are, however, affected by the separation time of two flares. When the separation interval 
is longer than the electron escape time (Stdyn)) the light curves can almost be regarded as a simple 
sequence of two separate single injections. But when the separation interval is shorter than the 
electron escape time, multiwavelength light curves become rather complicated. The main physical 
reason behind this complication is the dynamic accumulation of both photons and relativistic 
electrons. First of all, 1 eV and 1 keV emissions are from synchrotron and others are from SSC. 
All the SSC emissions have second peak higher than the first one, this is due to the increase both 
in soft photon energy density and in the number of electrons which are not completely cooled yet 
when the second injection occurs. This accumulation of electrons also accounts for the increase 
in 1 eV synchrotron flux. The 1 keV emissions, however, show a lower flux in the second peak, 
this is because the relativistic electrons from the second injection are subject to a much stronger 
cooling due to the enhanced photon energy density from the first injection. In addition to the 
flux differences, there are obvious delays in reaching the peak fluxes for different wavelengths with 
respect to the synchrotron and SSC peaks, though 1 keV and 1 TeV fluxes track each other rather 
well. 

As demonstrated in these figures, the slow response and relatively small amplitude variations 
at the photon energies other than the synchrotron and SSC peaks argue against the usual belief of 
closely correlated variations in multiwavelength observations. Only the emissions from the tail of 
the electron energy distribution can be reliably used as diagnostics for separate injections. Still, 
extra caution is obviously needed when relating the energy contained in nonthermal particles 
versus the observed fluxes. 

4. SUMMARY AND DISCUSSIONS 

Using a homogeneous synchrotron self-Compton model, we have calculated the time evolution 
of emission spectra and electron energy distributions when nonthermal electrons are uniformly 
injected into a relativistically moving plasma blob with constant velocity. We have found that: 

(1) When the luminosities of the synchrotron and SSC peaks are comparable, the electron 
cooling by inverse Compton scattering is not negligible and the system is inherently nonlinear and 
dynamic. One has to solve the time-dependent, coupled electron and photon kinetic-equations 
self-consistently. Furthermore, since observations are most sensitive to the peak fluxes of 
synchrotron and SSC components, accurate treatments of synchrotron emissivity due to the end 
point effects and inverse Compton scattering in the KN regime are quite essential. 

(2) When the cooling time at the maximum particle energy is longer than the injection 
timescale (> R/c), the light curve of emissions corresponding to the tail of the electron distribution 
can have short, large amplitude variations but emissions at other wavelengths show considerably 
longer and smaller amplitude changes. Additionally, spectral evolution is also rather slow. All 
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these features are simply caused by the long cooling time of electrons. 

(3) When cooling time at the maximum particle energy is shorter than the dynamic timescale, 
strong spectral evolutions are observed for both synchrotron and SSC components and short 
duration flares are obtained in most energy bands. 

(4) Generally, the prompt TeV spectrum is curved due to the KN effect and the fact that 
TeV-production electrons are usually in a cooled distribution. This consideration does not take 
into account the possible infrared background attenuation of the TeVs, which might cause further 
curvature in the TeV spectrum. On the other hand, most current TeV observations require an 
accumulation time probably much longer than the dynamic timescale of the blob, so that it 
might still be possible to obtain a quasi power-law TeV spectrum by averaging over an evolving 
spectrum. Further study is needed to address this issue. 

(5) We recommend plotting the light curves in a fashion that is logarithmic flux versus linear 
time. The goal is to look for exponential decays at specific energy bands, which might give direct 
measurements of the size of the system, as indicated in Figures ^-11. 



(6) One has to be cautious about the common belief that light curves in different energy bands 
should track each other. The electrons responsible for producing specific energy photons might 
have quite different lifetimes, especially when multiple and closely spaced injections are involved. 
This complication also applies to the leading/lagging analysis for different photon energy bands. 

(7) When high time-resolution spectroscopy is available both in keV and TeV bands, one 
should be able to prove whether TeV production is via SSC process by comparing the rates of 



spectral softening as done in Figure 15 



The primary purpose of this paper is to investigate the radiative signatures in a purely 
cooling and dynamic system, thus providing a bridge between observations and the detailed but 
largely unknown physics of particle energization processes. Since we did not address the particle 
acceleration problem here, in this sense, some of the conclusions drawn above are certainly subject 
to revisions as our understanding of energy flow in AGNs improves. 

In conclusion, we have found that solving time-dependent, coupled electron and photon 
kinetic-equations provides an easy and efficient way of comparing multiwavelength, time-dependent 
observations with some simplified SSC models. It has the advantage of naturally combining the 
spectral and temporal evolutions in a dynamic system, which is very useful when more and more 
high quality data become available. 
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Photon Energy (MeV) 

Fig. 1. — Time evolution of electron energy distribution (upper panel) and the corresponding photon 
spectra (lower panel) using parameters from the unique solution discussed in § |3.1.1 . Electrons are 
injected at a constant rate lasting 2tdyn in the comoving frame. There are 20 curves in each panel, 
which start at t = O.Stdyn and end at t = lOtdyn with a time interval of O.Stdvn between them. 



The injection process is shown by the dashed curves moving up in n{E), along with the increasing 
photon fluxes. After the injection stops, electrons are continuously cooled and the photon spectrum 
softens. These parameters give a comparable peak fluxes for synchrotron and SSC components. 
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Fig. 2. — Multiwavelength light curves in observer's frame (T> is assumed to be 10) using parameters 
from the unique solution (cf. Figure |l|). Fluxes at the synchrotron and SSC peaks show fast time 
variability with large amplitudes, but fluxes at other wavelengths have a very long plateau with 
very small amplitude variation. 
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Fig. 3. — Shown is |7ssc/7syn| as a function of electron energy at different times (from 0.5 — lOtdyn) 
using parameters from the unique solution (of. Figure |^). The horizontal solid line at the ratio 
being 1 is plotted to guide the comparison. The SSC process becomes important as soon as the 
photon energy is built up and, in fact, is more important than synchrotron cooling for most of the 
electron energies. 
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Fig. 4. — Shown is the photon spectral evolution from — 4idyn as the magnetic field is varied from 
the unique solution value Bg ■ The ratio r/ is unity when B = Bg (cf . Figure ||) but varies nearly as 
(Bg/B)^ as shown in the middle and lower panels here. 
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Fig. 5. — The photon spectra vF,^ for 6 different total-injection-energy £, ranging from 10 ergs 
to 10^^ ergs with a factor of 10 increment in each case. All 6 spectra are taken at the end of 
electron injection {t = 2tdyn) and their fluxes are divided by £^/10^^ so that if the synchrotron flux 
was exactly proportional to £, they would have had the same heights. As £ increases, the electron 
cooling changes from synchrotron dominated (when £ < 10^^) to synchrotron sclf-Compton (SSC) 
dominated. Note that the increase in SSC component is proportional to £^ for £ < 10^^. When 
SSC cooling is very strong, electrons cool so quickly that the synchrotron flux at 2t(jyn is no longer 
scaled as £ anymore as shown in the 10^^ (dashed) and 10^^ (dotted) cases. Also in these large £ 
cases, the efficient cooling makes the synchrotron peak rather broad. 
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Fig. 6. — Time evolution of electron energy distribution (upper panel) and the corresponding 
photon spectra (lower panel) for £ = 10^"^ ergs. Electrons are injected at a constant rate lasting 
2tdyn in the comoving frame, with 7inin = 10, 7max = 10^, and index p = 2. There are 20 curves in 
each panel, which start at f = O.Stdyn and end at t = lO^dyn with a time interval of 0.5tdyn between 
them. The injection process is shown by the dashed curves moving up in n(£'), along with the 
increasing photon fluxes. After the injection stops, electrons arc continuously cooled and photon 
spectrum softens (solid curves). With this injection energy, the SSC peak is considerably lower 
than the synchrotron peak. 
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Fig. 7. — Same as Figure |6|, except tliat £ = 5 x W^"^ ergs. Now tlie SSC component lias become 
comparable to the synchrotron component, and the whole system understandably evolves on a 
faster timescale. 
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Fig. 8. — Same as Figure ^, except that £ = 10^^ ergs. The SSC component dominates over the 
synchrotron component, and the buildup of synchrotron photon energy-density is so quick that 
electron cooling is very efficient. Towards the end of simulation (~ lOtdyn), multiple Compton 
scattering features are evident. 
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Fig. 9. — Multiwavelength light curves in observer's frame is assumed to be 10) with £ = 10 
ergs (cf. Figure ^). Solid and dashed curves are for synchrotron and SSC components, respectively. 
The exponential decay depicted by the keV and TeV fluxes allows a direct estimate of the size of 
the emission cloud. 
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Fig. 10. — Same as Figure ^ except that £" = 5 x 10^'^ ergs. (Also cf. Figure 



- 28 - 




-29- 




Oe+OO 1e-14 Oe+00 5e-18 Oe+00 2e-18 
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Fig. 12. The time evolution of the correlation between the photon index and the flux (ergs 

cm~^ s^^) at 2 kcV, 100 keV, and 1 GeV, respectively. Arrows indicate the direction of the time 
evolution. The total injected energy is £ = 10^^ ergs. A large spectral evolution is seen at 100 
keV, where synchrotron and SSC components mix. Spectral evolution at keV and GeV bands are 
relatively moderate and clockwise. 
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Fig. 15. — Evolution of peak photon energies of both synchrotron and SSC components. The total 
injected energy is <S = 5 x 10^^ ergs (cf. Figure |^. The upper panel shows their actual values 
and the lower panel shows the normalized values after dividing each peak energy by the value 
at t = O.Stdyn for synchrotron and SSC components, respectively. The synchrotron peak energy 
(oc 7^) decreases faster than SSC peak energy (oc 7) initially, because SSC process is still in the 
Klein-Nishina regime. 
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Fig. 16. — Light curves of observed spectra at 1 kcV when power-law electrons are injected. Solid 
curve is calculated when Q{'y) = <3o7^^ for < t < 2tdyn, and Q(7) = for t > 2tdyn- Dotted 
curve is obtained when Q(7) = Qo7~^i/idyn for < t < 2tdyn, and Q(7) =0 for t > 2tdyn- The 
upper panel shows the time profile of the injection Q in arbitrary units. 
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Fig. 17. — Expected light curves at different wavelengths (lower panel) when power-law electrons 
are injected according to the profile in the upper panel. The time is measured in the observer's 
frame. In the comoving frame, two injections are separated by Stdyn, longer than the electron 
escape timescale, which is chosen as Stdyn- The same amount of energy, 5 x ergs, is injected 
in each flare. The dotted, dashed, dash-dotted, dash-dot-dot-dotted, and solid represent fluxes at 
1 eV, 1 keV, 1 MeV, 1 GeV, and 1 TeV, respectively. The two flares can be regarded as a simple 
sequence of two unrelated injections. 
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Fig. 18. — Same as Figure 17, except that the separation of two flares is shorter (2tdyn) than the 
electron escape timescale. The second flare is now strongly affected by the residual effects from the 
first electron injection. 



